Electronic structure of self-assembled InAs/InP quantum dots: A Comparison with 

self-assembled InAs/GaAs quantum dots 
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We investigate the electronic structure of the InAs/InP quantum dots using an atomistic pseu- 
dopotential method and compare them to those of the InAs/GaAs QDs. We show that even though 
the InAs/InP and InAs/GaAs dots have the same dot material, their electronic structure differ 
significantly in certain aspects, especially for holes: (i) The hole levels have a much larger energy 
spacing in the InAs/InP dots than in the InAs/GaAs dots of corresponding size, (ii) Furthermore, 
in contrast with the InAs/GaAs dots, where the sizeable hole p, d intra-shell level splitting smashes 
the energy level shell structure, the InAs/InP QDs have a well defined energy level shell structure 
with small p, d level splitting, for holes, (iii) The fundamental exciton energies of the InAs/InP 
dots are calculated to be around 0.8 eV (~ 1.55 /im), about 200 meV lower than those of typical 
InAs/GaAs QDs, mainly due to the smaller lattice mismatch in the InAs/InP dots, (iii) The widths 
of the exciton P shell and D shell are much narrower in the InAs/InP dots than in the InAs/GaAs 
dots, (iv) The InAs/GaAs and InAs/InP dots have a reversed light polarization anisotropy along 
the [100] and [110] directions. 

PACS numbers: 68.65.Hb, 73.22.-f, 78.67.Hc 
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I. INTRODUCTION 

Self-assembled semiconductor quantum dots (QDs) 
have attracted a large interest because the discrete 
and isolated energy levels due to the 3D confinement 
can be utilized in high-efficiency and low-threshold cur- 
rent lasersjiSi&i single photon emitters 5 -! 6 - and in quan- 
tum computing applications (qubits), 7,8 1.3 /.tm opti- 
cal devices have been achieved using InAs/GaAs quan- 
tum dots^i 3 ^ however, it is difficult to push further 
the InAs/GaAs QD devices to work at the more de- 
sirable ^1.55/im telecommunication wavelength because 
the large compressive strain (7%) in the InAs/GaAs dots 
enlarges the conduction-valence energy gap of InAs too 
much for this purpose. Intuitively, a InAs dot embed- 
ded in a less lattice mismatched host material, e.g. InP, 
should solve the problem. Indeed, recently there have 
been reports on 1.55 /mi InAs/InP QDs lasers with a 
high-gain and a low-threshold current which opens 
up the possibility to integrate quantum dot material into 
an optical cavity^ 1 - and hence to reliably fabricate a sin- 
gle photon source in the telecommunications wavelength 
ranged 

Despite their importance, there are only few stud- 
ies on the InAs/InP system, both theoreticall y 12 ' 13 and 
experimentall y 14 ' 15 ' 16 ' 17 ' 18 ' 19 compared to the well stud- 
ied InAs/GaAs QDs. Although the InAs/InP dots and 
the InAs/GaAs dots have the same dot material, they 
differ in three aspects, (i) InAs/InP has a much smaller 



lattice mismatch (3%) than InAs/GaAs (7%). (ii) The 
InAs/InP dots have a less confining potential for elec- 
trons, but a stronger confinement for holes than the 
InAs/GaAs dots, (iii) The InAs/InP dots share the same 
cation (In), while the InAs/GaAs QDs share the same 
anion (As) at the interface. These differences may lead 
to different electronic and optical properties of the two 
systems. The continuum theories, such as those based 
on the effective mass approximation (EMA) and multi- 
bands kp method, may in principle capture the first two 
differences, however, to account for the third aspect, one 
needs atomistic theories such as the empirical pseudopo- 
tential method s 20 ' 21 or the tight-binding methodsj 12 ' 22 
The empirical pseudopotential methods have been suc- 
cessfully applied to various systems. 21 ' 23 In this paper, 
we perform a comparative study on the InAs/InP QDs 
and InAs/GaAs QDs using an atomistic pseudopoten- 
tial method. We find that there are significant differ- 
ences in the electronic structure between the two sys- 
tems, including the single-particle energy levels and opti- 
cal properties. These differences, which have not yet been 
paid enough attention to, could be revealed in the fu- 
ture high-resolution optical spectroscopy 24 ^ and charg- 
ing experiments 



26,27 
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The rest of the paper is organized as follows. In Sec. 
III1 we introduce briefly the atomistic pseudopotential 
method used in the calculations. In Sec. IIII[ we compare 
the strain profiles and the strain modified band-offsets 
for the InAs/GaAs and InAs/InP QDs. In Sec. [TV] we 
compare the pseudopotential calculated electronic struc- 
ture as well as the wavefunctions of the two dots. We 
compare the excitonic transitions of the two dots in Sec. 
IVl and conclude in Sec. IVII 
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II. METHODS 



The geometries of the QDs studied here are lens- 
shaped InAs dots embedded in the host materials (GaAs 
or InP) matrices containing 60 x 60 x 60 8-atom unit 
cells. The dots are assumed to grow along the [001] di- 
rection. We performed the calculations on dots with base 
diameters D= 20, and 25 nm, and for each base diameter, 
we vary the dot heights h from 2.5 nm to 5.5 nm, as the 
dot height is relatively easy to control in experiments. 28 
We first relax the dot+matrix system by minimizing the 
strain energy as a function of the coordinate {R„ !Q ,} of 
atom a at site n for all atoms, using valence force field 
(VFF) methodsi 29 i 30 Once we have the atom positions, 
we obtain the energy levels and wavefunctions by solving 
the single-particle Schrodinger equation, 



\V 2 + V ps {v) 



ipi(r) = e f Vi(r) 



(1) 



where Vp S (r) = V so + J2 n J2 a v a( r ~ R-n,a) is the su- 
perposition of local screened atomic pseudopotentials 
v a (r), and the total (non-local) spin-orbit (SO) poten- 
tial Vso- The screened pseudopotentials^ are fitted to 
the physical important properties of the materials, in- 
cluding the band energies at high-symmetry points, ef- 
fective masses, strained band offsets, hydrostatic and bi- 
axial deformation potentials of individual band edges. 
The pseudopotentials of InAs/InP dots are given in Ap- 
pendixO whereas the potentials for InAs/GaAs QDs are 
taken from Ref. [U. The Schrodinger Eq. |T]) is solved 
by expanding the wavefunctions ipi as linear combina- 
tions of bulk bands (LCBB)2£ {(j) m ^^(k)} of band in- 
dex m and wave vector k of material A (= InAs, GaAs, 
InP), strained uniformly to strain *e. We use m = 2, 3, 4 
for the hole states, and m = 5 for electron states on a 
6x6x16 k-mesh. We use *e = for the matrix mate- 
rial, and an average *e value from VFF for the strained 
dot material (InAs). It has been shown that the en- 
ergy levels changes in InAs/GaAs QDs due to the piezo- 
electric effects are quite small^i Because the lattice mis- 
match in the InAs/InP QDs is only half of that of the 
InAs/GaAs QDs, we expect that the piezo-effect should 
be even smaller in the InAs/InP dots, and therefore is 
ignored in the present calculations. 

The exciton energies are calculated using the configu- 
ration interaction (CI) method following Ref. 



III. STRAIN PROFILES AND STRAIN 
MODIFIED BANDOFFSETS 
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FIG. 1: Strain profiles and strain modified band-offsets for 
lens-shaped InAs/GaAs QD (D=25 nm, /i=3.5 nm). The 
strain profiles are shown in (a) along the [100] direction and 
in (b) along the [001] direction. The strain modified band- 
offsets are shown in (c) along the [100] direction and in (d) 
along the [001] direction for the CBM, HH, LH and SO bands. 
The reference energy is chosen to be the VBM of GaAs. 



03 
CO 




^Hydrostatic strairj 

IV. ,\ 



1 1 1 
" InP 

'(b) 


! lnP ' 








-15-10 -5 5 10 15 
Along [100] (nm) 



-15-10 -5 5 10 15 
Along [001] (nm) 



We first compare the strain profiles in the InAs/GaAs 
and the InAs/InP QDs. Figure Q] (a), (b) depict the 
hydrostatic and biaxial strain along the [100] direction 
and the [001] direction respectively for the InAs/GaAs 
dot, whereas Fig. [3] (a), (b) depict the strain pro- 
files for the InAs/InP quantum dots. The strain is cal- 
culated for the lens-shaped QDs with D=25nm, and 



FIG. 2: Strain profiles and strain modified band-offsets for 
lens-shaped InAs/InP QD (D—25 nm, h=3.5 nm). The strain 
profiles are shown in (a) along the [100] direction and in (b) 
along the [001] direction. The strain modified band-offsets 
are shown in (c) along the [100] direction and in (d) along 
the [001] direction for the CBM, HH, LH and SO bands. The 
reference energy is chosen to be the VBM of InP. 



3 



/i=3.5nm. The hydrostatic(isotropic) strain is denned as 
I = Tr(e) = e xx +£ yy +£ zz , reflecting the relative change of 
volume i.e., I ~ AV/V, whereas the biaxial strain is de- 
fined as B =y/{e xx - e vv ) 2 + {e zz - e xx ) 2 + {e yy - e zz ) 2 . 
As we can see, the hydrostatic strains of the InAs/InP 
and InAs/GaAs dots have very similar features: both 
are almost constants inside the dots and decay rapidly to 
zero in the matrices. The hydrostatic strain is negative 
in both dots, suggesting that the InAs dot is compressed. 
However, the hydrostatic strain in the InAs/InP dot is 
about half of that in the InAs/GaAs dot. The biaxial 
strain of both QDs also has similar features, decaying 
fast in the [100] and [010] directions outside the dots, it 
has long tails along the [001] direction. 

The hydrostatic strain and the biaxial strain will shift 
the conduction band minimum (CBM) and the valence 
band maximum (VBM), and therefore modify the band- 
offsets of the QD. The biaxial strain will further split 
the heavy-hole (HH) and light-hole (LH) bands. We 
analyze the strain-modified band-offsets via Pikus-Bir 
model ) 33 ' 34 i 35 ' 36 using the local strain input from VFF 
calculations. This model, however, serves only as an 
illustration of the strain effect and is not used in our 
actual calculation of the single particle energy levels. 
The strain modified band-offsets are illustrated in Fig. 
Q] (c),(d) for InAs/GaAs quantum dot, and in Fig. [2] 
(c),(d) for InAs/InP quantum dot along the [100] direc- 
tion and [001] direction respectively. The VBM of the 
unstrained host materials is set to be zero as the ref- 
erence energy. The band-offsets in the InAs/GaAs dot 
are 480 meV for the electron and 280 meV for the hole, 
compared with 1050 meV for the electron and 46 meV 
for the hole in the unstrained system^ The confinement 
potential for electrons is stronger than that for holes. In 
contrast, for the unstrained InAs/InP system, the band- 
offset^ is 580 meV for the electron and 420 meV for 
the hole, which change to 320 meV and 530 meV respec- 
tively in the InAs/InP QDs, due to the strain effects. 
The confinement potential for holes is stronger than that 
for electrons^ Therefore, the confinement potentials in 
the InAs/GaAs and the InAs/InP QDs are very differ- 
ent. While the confinement for electrons is weaker in 
the InAs/InP QDs, the confinement for holes is signifi- 
cantly stronger. How the different confinement potentials 
lead to different electronic structure in the two dot sys- 
tems will be discussed in Sec. IIV1 It was pointed out 
in Ref. HH that in the tall InAs/GaAs QD, the biax- 
ial strain might develop hole traps that localize holes at 
the interface of the QDs. This is unlikely to happen in 
the InAs/InP QDs for two reasons, (i) The strain in the 
InAs/InP QDs is much smaller than in the InAs/GaAs 
QDs. (ii) The band-offset for holes in the InAs/InP QDs 
is much larger than in the InAs/GaAs QDs. 

The biaxial strain splits the HH and LH bands in ad- 
dition to shifting the VBM. HH is higher in energy than 
the LH band, i.e., E^h > Eih inside both dots. The HH- 
LH splitting is about 180 meV for the InAs/GaAs QD 
and 120 meV for the InAs/InP QD. Outside the QD, the 



heavy- light-hole splitting changes sign, i.e., Ehh < Eih- 35 
It is also interesting to note that the SO band inside the 
InAs/GaAs quantum dot is lower in energy than out- 
side the QDs, while the opposite is true for the InAs/InP 
quantum dot, due to the large band-offsets between InAs 
and InP for holes. 



IV. SINGLE-PARTICLE ENERGY LEVELS AND 
WAVEFUNCTIONS 

In this section, we compare the pseudopotential calcu- 
lated single-particle energy levels as well as the wavefunc- 
tions of the InAs/InP QDs to those of the InAs/GaAs 
QDs. The energy levels are compared in three scales: (i) 
The energy difference between the lowest electron state 
and highest hole state, which largely determines the exci- 
ton energies and is in the order of 1 eV, (ii) The intraband 
energy spacing which is in the order of a few tens meV. 
(iii) The intraband p level splitting due to the G2 V atom- 
istic symmetry of the QDs, which is in the order of a few 
meV. The single-particle energy levels of the InAs/InP 
dots are summarized in Table. Q] whereas the results for 
the InAs/GaAs dots can be found in Table I of Ref. [H. 
The results of the InAs/GaAs QDs are very similar to 
what was obtained in Ref. l2ll 



A. Confined states and wavefunctions 

Figure [3fa) depicts the energy levels of all confined 
electron states and 6 highest confined hole states of the 
InAs/InP QDs for various sizes. As a reference, we show 
in Fig. [UJb) the lowest 6 confined electron levels and 
the highest 6 confined hole levels of a InAs/GaAs QD of 
D=20 nm, h=2.5 nm. The zero energy is chosen to be 
the VBM of the host materials. The confined electron 
(hole) states are defined to be the states whose ener- 
gies are lower (higher) than the CBM (VBM) of the host 
materials. We also show the electron and hole envelope 
wavefunctions of the InAs/GaAs and InAs/InP QDs for 
two dot geometries: a flat dot (D=25 nm, h=2.5 nm) 
in Fig. g] (a) (b) and a tall dot (D=25 nm, h=5.5 nm) 
in Fig. Q] (c) (d). We show the lowest 6 electron and 
the highest 6 confined hole states for each dot. In all 
cases, the isosurface is chosen to enclose 50% of the total 
charge. The number on the bottom of each small panel 
gives the percentage of the density localized inside the 
QD. 



1. Confined electron states 

As we can see from Fig. [3j the InAs/InP QDs 
have fewer confined electron states compared to the 
InAs/GaAs QDs due to the smaller conduction band off- 
set. In the smaller InAs/InP QD with D=20 nm, h=2.5 
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FIG. 3: Single-particle energy levels for (a) the lens-shaped 
InAs/InP QDs of different sizes, and (b) the lens-shaped 
InAs/GaAs QD (D=20 nm, h=2.5 nm). We show all con- 
fined electron states and the six highest confined hole states. 
The zero energies are chosen to be the VBM of the host ma- 
terials. 



nm, only one electron state (eo) is confined. When we in- 
crease the dot base and height, more states are confined. 
For the dot with D=25 nm, h=2.5 nm, eo, &i and e2 are 
confined, whereas in the dot with D=25 nm, /i=3.5 nm, 
the e4 and e$ levels are also confined. The confined elec- 
trons show well defined s,p, d energy level shell structure 
in both the InAs/InP and the InAs/GaAs dots. 



For the InAs/GaAs QDs, the electron wavefunctions 
of the flat dot [See Fig. H (a)] and the tall dot [See Fig. 
H] (c)] have very similar shapes. The lowest state has s- 
like wavefunction, followed by 2 p-like states. The lower 
energy p orbital (ei) has peaks along the [110] direction, 
whereas the higher energy p orbitals (e2) has peaks along 
the [110] direction^ Considering the next three levels, 
close in energy, the first two orbitals (e^, e$) that have 
similar shapes are d orbitals, whereas the third orbital 
(es) having a peak in the dot center, is the 2s orbital. 
The wavefunctions of the InAs/InP dots have very similar 
shapes to the corresponding ones of the InAs/GaAs dots, 
except that the wavefunctions of the InAs/InP dots are 
larger, indicating that the electrons are less confined in 
the InAs/InP dots. In the flat InAs/InP QD, only about 
53% of eo state is confined in the dot compared to 78% 
in the InAs/GaAs QD of the same size. Even in the tall 
InAs/InP QD, less than 82% of the eo state is confined 
compared to 92% in the tall InAs/GaAs QD. 



2. Confined hole states 

The valence band offset for the InAs /InP QDs is about 
530 meV, which is almost twice as much as in the 
InAs/GaAs QDs. Therefore, more hole states are con- 
fined in the InAs/InP dots than in the InAs/GaAs dots. 
Unlike the InAs/GaAs dot, where the energy level shell 
structure for holes are not so obvious^ the InAs/InP 
dots have well defined hole s, p, d energy level shell struc- 
ture, similar to those of electrons. 

The shape of the hole wavefunctions in the InAs/GaAs 
dot is more complicate than that of the electron wave- 
functions. But nevertheless in the flat dot [See Fig. 0(a)] 
they can still be recognized as s, p, d, 2s orbitals although 
with some mixed characters^ In the tall InAs/GaAs dot, 
the holes are strongly localized on the interface of the dot 
due to the strain effects^ and do not have clear s, p, d 
characters any more [See Fig. H](c)]. On the other hand, 
the hole wavefunctions of the InAs/InP dots are quite 
different from those of the InAs/GaAs dots. They are 
quite similar to the electron wavefunctions, for both the 
flat and tall dots, except that the two p orbitals switch 
order in energy, i.e., the first hole p orbital has peaks 
along the [110] direction, whereas the second hole p or- 
bital has peaks along the [110] direction. The rotation of 
the p orbitals has also been noticed by Sheng et al.^2 and 
was attributed to the piezoelectric effects. However, in 
our calculations, the piezo-effect is ignored, and therefore 
it cannot be the reason for the p orbital rotation. The 
rotation of the wavefunctions can be explored experimen- 
tally via magnetocapacitance spectroscopy. 40 Unlike in 
the InAs/GaAs dots, no hole localization has been found 
in the InAs/InP dots, because of their smaller strain and 
the large confinement potential for holes, as is discussed 
in Sec. [ml 

The holes are strongly confined in the InAs/InP QDs. 
Even in a small InAs/InP QD, more than 80% of charge 
density is localized in the QD for the ho and hi 2 states. 
As a result, in a small InAs/InP QD, electrons could es- 
cape from the dot much more easily than holes, resulting 
in a positively charged system. Therefore the InAs/InP 
dot might be a good candidate for a memory device via 
hole storage in the doti 14 i 41 ' 42 



3. Single-particle electron-hole energy gap 

The single-particle electron-hole energy gap Ae e- h — 
eo — ho is summarized in Table HI For the InAs /InP QDs 
with D=20 nm, and h=2.5 nm, Ae ej ^ ~ 917 meV, which 
is about 180 meV less than that of the InAs/GaAs dot of 
the same size. The electron energy levels move down in 
energy, whereas the hole energy levels move up in energy, 
with the increasing of the base and height of the dot, as 
shown in Fig. [3] (a). As a result, Ae Ct h decreases with the 
increasing of the dot size, due to the reduced confinement. 
For the dot of D=20 nm, Ae e> /, changes from 918 meV 
at h~2.5 nm, to 798 meV at h—5.5 nm, with a 120 meV 
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FIG. 4: (Color online) Top view of the squared wavefunctions of the confined electron and hole states for the InAs/GaAs and 
InAs/InP QDs. We show in (a), (b) the results of flat dots with D=25 nm, h=2.5 nm, and in (c) and (d) the results of tall 
dots with D=25 nm, h—5.5 nm. The isosurface to chosen to enclose 50% of the density. The number on the bottom of each 
small panel gives the percentage of the density localized inside the QD. 



reduction. A similar Ae e ^ reduction has been found for 
the InAs/InP dots of base D=25 nm when the dot height 
increases from 2.5 nm to 5.5 nm. The electron-hole en- 
ergy gap reduction results in a corresponding redshift of 
the exciton emission lines (Sec. |V|) . 



B. Intraband energy spacings 

The intraband energy spacing can be used to charac- 
terize the confinement effects in the QDs. We summarize 
the s-p energy spacing S sp = e± — eo (or ho — hi), and 
p-d energy level spacing 8 p d — &z — e.2 (or h% — ft.3) of 
the InAs/InP dots in Table []] whereas the results for the 
InAs/GaAs dots can be found in Table I of Ref. [H. To 
further see the trend of how the energy spacing changes 
with the dot size, we plot 8 sp of both InAs/GaAs and 
InAs/InP QDs, as a function of the dot height in Fig. [5] 
(a), (b) for electrons and holes respectively. 

(a) electrons: The electron energy spacing of the 
InAs/InP dots is in the range of 50 to 70 meV, slightly 



smaller than the energy spacing (50 to 80 meV) of the 
InAs/GaAs dots, due to the weaker confining potential 
for electrons in the InAs/InP dots. The s-p energy spac- 
ing S sp and p-d energy spacing S p d are nearly equal, in 
rough agreement with the EMA with harmonic confine- 
ments. 

Intuitively, S sp should decrease monotonically by in- 
creasing the dot height. This trend is followed by the 
InAs/GaAs dots. Surprisingly, the electron 5 sp of the 
InAs/InP dots increases with the increasing of the dot 
height, against the naive expectation. In this case, both 
s, p levels move down in energy with the increasing of 
the dot size, but the s level moves down faster than the 
p levels ^ leading to a larger 5 sp . 

(b) holes: The hole energy spacing in the InAs/InP 
dots ranges from 20 to 40 meV, significantly larger than 
that (< 20 meV) of the InAs/GaAs dots of the same 
size. This is because the hole confining potential in 
the InAs/InP dots is about 340 meV larger than in the 
InAs/GaAs dots. Nevertheless, in all the InAs/InP QDs 
we have studied, the electron energy spacing is still about 
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TABLE I: Summary of the pseudopotential-calculated single- 
particle level spacing (in meV) of the InAs/InP quantum dots 
of different base sizes and heights (in nm). Unconfined states 
are leave in blank in the table. 
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FIG. 6: Comparison of the p level splitting S pp between 
the InAs/InP and InAs/GaAs QDs for (a) electrons, and (b) 
holes. 
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C. Intraband p level splitting 

In a continuum model, ignoring the underlying atom- 
istic structure, a cylindrical QD has the Coo V symmetry, 
leading to degenerate p, d levels. In contrast, the atom- 
istic theories maintain the real symmetry of the QDs. 
For example, the lens-shaped QDs made of zinc-blende 
III-V semiconductors are of Ci v symmetry,— where the 
[110] and [110] directions are non-equivalent, resulting in 
split p levels and d levels. The values of the p-p split- 
ting Spp — e2 — ei (or hi — /12) are summarized in Table 
U for the InAs/InP QDs. We also depict the p-p split- 
tings as functions of dot height for both InAs/InP and 
InAs/GaAs QDs in Fig. H 



FIG. 5: Comparison of the intraband energy spacing 8 3p be- 
tween the InAs/InP and InAs/GaAs QDs for (a) electrons, 
and (b) holes. 



twice larger than that of the holes, because electrons have 
a much lighter effective mass than holes. For holes, 5 p d 
is much smaller than 5 sp , deviating from the harmonic 
potential approximation. As shown in Fig. [5] (b), the 



p level splitting for electrons 



The electron p-p splittings are shown in Fig. [5] (a) as 
functions of the dot height. The electron 5 PP is ~ 1 - 
2 meV for flat InAs/InP dot, and increases to about 3 
- 5 meV for the tall dots. The dots with the smaller 
base (D=20 nm) show a much larger p-p splitting. The 
electron p-p splitting for the InAs /InP dots is about 3%- 
5% of the electron s-p energy spacing S sp . The electron 
p-p splittings for the InAs/GaAs dots are in the range 
of 2 - 4 meV, and 3 % - 4% of S sv , and show a weak 



energy spacing 5 sp of holes decrease monotonically with dependence on the dot height, 
the increasing of the dot height for both InAs / GaAs and 
InAs/InP QDs. No anomaly is found. Notice that for 
very tall QDs, S sp becomes very small for the InAs/GaAs 
QDs, due to the hole localization on the interface of 
the QDs [see Fig. @] (c)], which is not the case for the 
InAs/InP QDs. 



2. p level splitting for holes 

The hole p-p splittings are shown in Fig. [S] (b). It has 
been shown in previous studies ; 21 ! 38 ' 44 that the hole p-p 
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FIG. 7: Comparison of the primary exciton energies vs. dot 
height between the InAs/InP and InAs/GaAs QDs. 



splittings can be larger than 10 meV in the InAs/GaAs 



sp • 



dots, which is about 70% of the s-p energy spacing S, 
For very tall InAs/GaAs dots, in which the holes lo- 
calize on the interface, 8 PP could even be much larger 
than S sp . The large p-p splitting leads to very differ- 
ent electronic and optical properties of the InAs/GaAs 
QDs than those predicted by continuum theories, e.g. 
the nontrivial charging pattern that breaks Hund's rule 
and the Aufbau principle for holes ; 38 ' 44 that has been 
recently confirmed experimentall y 27 ' 45 Surprisingly, the 
calculated hole p-p splitting of the InAs /InP dot is much 
smaller than that of the InAs/GaAs dots, even though 
they have the same dot materials. For the flat InAs/InP 
dots, the splitting is only about 1 ~ 2 meV and for the 
tall dot with a small base (D =20 nm) the splitting can 
be as large as 4.2 meV. Nevertheless, the hole S pp is less 
than 10 % of 5 sp for the InAs/InP dots. Therefore the 
multi-hole phase diagram and the charging patterns in 
the InAs/InP dots are expected to be very different from 
those in the InAs/GaAs dots, which can be examined by 
the hole charging experiments^ This is one of the most 
important results of the present work, which can not be 
obtained from the continuum theories. For example, a 
k-p theory predicted small (~ 1 meV) hole p-p splitting 
for both types of dotsj 12 ' 46 ' 47 The difference in the hole 
splitting for the two types of dots might comes from their 
different strain profiles, the band offset, or the interface 
effects (common cation vs. common anion). We leave 
this for future investigations. 



V. EXCITONS 

Figure [7] depicts the fundamental exciton energies vs. 
dot height for the InAs/InP and InAs/GaAs dots. The 
exciton energies of the InAs/GaAs dots are about 200 
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FIG. 8: (a) The exciton transition strengths for the InAs/InP 
QD of D=25 nm, h=3.5 nm. The upper panel is for the 
transitions polarized in the [110] direction, while the lower 
panel is for transitions polarized in the [110] direction. X 3 
is the polarization anisotropy (Eq. [3} of the S exciton. (b) 
The exciton P shell width vs. dot height, (c) The S exciton 
polarization anisotropy X 3 vs. dot height. 



meV higher than those of the InAs/InP dots. The exci- 
ton energies decrease monotonically with increasing of 
the dot size. For the InAs/InP dot with L>=20 (25) 
nm, the exciton energy reduces from about 900 (870) 
meV to about 780 (750) meV, as the dot height increases 
from 2.5 nm to 5.5 nm. The exciton energy reduction is 
about 120 meV, compared to 70 meV for the InAs/GaAs 
dots. One of the most important motivations of study- 
ing the InAs/InP dots is the 1.55 /Ltm (~ 800 meV) 
emission for device applications. This wavelength can 
be easily achieved by the InAs/InP dots with a reason- 
able size, while it is challenging to be obtained using the 
InAs/GaAs dots, as shown in Fig. \7\ 

In the Hartree-Fock approximation, the fundamental 
exciton energy can be calculated as, 



E 



x 



j(eh) 



(2) 



where Ae e ^h — eo — ho is the single-particle electron- 
hole energy gap and j( eh > is the direct electron-hole 
Coulomb energy. The electron-hole exchange energy is 
much smaller than j( e ' 1 ) and therefore is ignored here. 
In typical InAs dots, Ae e ./j ^> j( eh K Therefore the ex- 
citon energy is largely determined by the single-particle 
electron-hole energy gap Ae e .h (See Sec: IIV A 31) . 

We further calculated the higher excitonic transitions 
for the InAs/InP dot with D—25 nm and /i=3.5 nm, 
for which the fundamental exciton wavelength is close 
to 1.55 /im. The results are shown in Fig. [5] (a), where 
the upper panel depicts the transition strength along the 
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[110] direction, while the lower panel shows the transition 
strength along the [110] direction. The transitions polar- 
ized along the [001] direction are about 4-5 orders of mag- 
nitude smaller than those parallel to the (001) plane, and 
are therefore not considered here. The excitonic transi- 
tions form several shells: the first shell coming from the 
cq to ho transitions is the exciton S shell. The transi- 
tions of ei.2 to hi t 2, form the exciton P shell, whereas 
the transitions of 634,5 to hz.i,h form the D shelh 25 i 48 
There is also a small transition peak at about 0.857 eV, 
which comes from the recombination of eo electron and 
/15 hole. This transition, however, is significantly weaker 
in the InAs/GaAs QDs. The exciton P shell and D shell 
are composed of a bunch of transitions with slightly dif- 
ferent transition energies. For the InAs/GaAs QDs, the 
width of the P shell is about 12 meV, whereas for the 
InAs/InP QD, the width is only about 2.5 meV. The D 
shell width is about 10 meV for the InAs/InP QDs, com- 
pared to about 22 meV in the InAs/GaAs dots. The P 
shell width with respect to the height of the QDs is pre- 
sented in Fig. M (b) for the dot with £>=25 ran. The P 
shell width of the InAs/GaAs increases to about 20 meV 
for h =5.5 nm. However, for InAs/InP dots, we found 
that the width is always about 2 ~ 3 meV. This feature 
reflects the fact that the p level splitting is small for both 
electrons and holes in the InAs/InP dots, fact that could 
be examined by the single-dot optical spectroscopy 24 ' 25 
We also calculated the light polarization anisotropy A, 
defined as the ratio of the transition intensities along the 
[110] and [110] direction, i.e.,21 

A = ^M. (3) 

I[U0] 

The results for the S shell excitons in the dot with 
base D=25 nm are presented in Fig. [5] (c) for both 
the InAs/InP and InAs/GaAs dots. For the InAs/GaAs 
dot, we found X s > 1, indicating that the intensity along 
the [110] direction is stronger than that along the [110] 
direction^ whereas, in the InAs/InP QD, A s < 1, indi- 
cating that the stronger intensity is along [110] direction. 
This feature could also be examined by the optical spec- 
troscopy. 



VI. CONCLUSION 

We have studied the electronic structure of the 
InAs/InP QDs using an atomistic pseudopotential 
method and compared them to those of the InAs/GaAs 
QDs. Our results show that even though the InAs/InP 
QDs and InAs/GaAs QDs have the same dot material, 

I 



V a (r ~ R na ) = v a (r - R na ) [1 + Sv na (ej] 



their electronic structure and optical properties differ sig- 
nificantly in certain aspects. These features, which may 
have important impacts for device applications, could be 
examined in future experiments. Some of the features 
can only be captured by atomistic theories and there- 
fore provide a unique opportunity to test the predictive 
capability of the different theoretical approaches. 

TABLE II: Band parameters obtained from the pseudopo- 
tential band structure and the target values of the fit. The 
"Target values" are conventional bulk parameters used in the 
literature (see Refl37l). AEvbo is the valence band offset 
relative to the bulk InAs VBM. Ao is the spin-orbit split- 
ting, m* are the effective masses at F, and a g , a v and b are 
the hydrostatic deformation potentials of the band gap, the 
valence band maximum, and the biaxial deformation poten- 
tial of the valence band, respectively. The predicted band 
structure critical points are compared with the existing ex- 
perimental data^. 



Parameters InAs InP 





PP 


Target 


PP 


Target 


Fit 


E g (eV) 


0.410 


0.410 


1.424 


1.424 


AEvbo (eV) 


-0.006 


0.000 


-0.440 


-0.420 


Ao (eV) 


0.390 


0.390 


0.109 


0.108 


* 


0.022 


0.024 


0.059 


0.080 


mJfc(OOl) 


0.387 


0.341 


0.444 


0.520 


m^(Hl) 


1.006 


0.917 


1.180 


0.950 


mjfc(001) 


0.027 


0.027 


0.085 


0.110 


m i \(Hl) 


0.026 


0.026 






ro: o (001) 


0.097 


0.085 


0.152 


0.21 


a g 


-6.44 


-6.6 


-6.93 


-6.0 


a v 


-1.01 


-1.0 


-0.68 


-0.6 


b 


-1.78 


-1.70 


-1.67 


-2.0 


Predictions 


T 7c (eV) 


4.55 


4.52 


5.31 


4.72 




-2.38 


-2.4 


-2.38 


-2.3 


X-7 V 


-2.37 


-2.4 


-2.25 


-2.2 


X% c 


2.28 




2.21 


2.38 


Xl c 


2.29 




2.61 




L§ v 


-1.14 


-0.90 


-0.92 


-1.23 


La,5v 


-0.87 


-0.90 


-0.80 


-1.12 


L(j C 


1.46 




2.15 


2.03 



APPENDIX A: EPM FOR INAS/INP 

As explained in Sec. II the crystal potential is written 
as a superposition of atomic potentials v a centered on the 
atomic positions. For each atomic potential we use for 
the screened pseudopotentials the expression proposed by 
Williamson et ali^i 



£ e^-^Wdql) [1 + Sv na (e)} , (Al) 
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TABLE III: Fitted pseudo-potential parameters for InAs/InP. 
A plane- wave cutoff of 5 Ryd is used. 



parameters As(In) P(In) In(As) In(P) 





56.8819 


0.1509 


853.4653 


5012.0545 


ai 


2.7023 


2.9215 


1.9724 


1.8556 




1.4894 


1.2190 


19.1236 


88.8570 


a 3 


0.5757 


0.3554 


0.5439 


0.7419 




0.00 


0.00 


1.6597 


1.6460 


Q.SO 


0.1315 


0.0140 


0.4056 


0.4800 



where u Q (|q|) has the functional form: 

2 

v a (\q\)=a 0a l a , (A2) 

a 2a e a3aq - 1 

and 

Sv na (e) = 7 Q • (e xx + e yy + e zz ). (A3) 

en are elements of the local strain tensor. The term 
5v na (e) plays a crucial role in describing the absolute hy- 
drostatic deformation potentials, in particular the vari- 
ation of the valence band edge and, separatly, the con- 
duction band edge under arbitrary strains. This allows 
us to describe the modification of the valence and con- 
duction band offsets when the systems are subjected to 
hydrostatic or biaxial deformation conditions such as in 
the case of epitaxial growth on a lattice-mismatched sub- 
strate. The parameters entering the previous equations 
have been determined by fitting a number of experimen- 
tally and theoretically (ab-initio) determined properties 
of bulk InP and InAs: the experimentally measured elec- 
tron and hole effective masses, band gaps (target values 
at 0° K), and spin-orbit splittings, hydrostatic deforma- 
tion potentials of the band gaps, band offsets, and LDA- 
predicted single band edge deformation potentials^. In 



the previous equation the term (3 has been introduced to 
represent the quasiparticle nonlocal self-energy effects. 
This kinetic energy scaling is needed to simultaneously 
fit bulk effective masses and band gaps. 

In Table. ITT1 we report the target values we have fit for 
the binary InAs, and InP, and the results of the fitting 
procedure. The target values correspond to the band 
parameters used in the literature^ at T = OK. A 5 Ry 
kinetic cutoff was used when generating the pseudopo- 
tentials. This cutoff has then been used in the QD cal- 
culations of this paper. From Table HT1 we see that the fit 
is satisfying. The corresponding parameters of the em- 
pirical pseudopotentials are given in Table IIIIl Although 
we fitted only a few band properties per material, we 
checked that the fit works also for the full band struc- 
ture. The predicted (not fitted) critical point energies 
are also reported in Table [TTJ Notice that we are using 
slightly different potentials for the In atoms in InP and 
InAs, to take into account the different charge redistri- 
bution occurring around the In atom when it is placed 
in a different environment. Considering only the nearest- 
neighbor environment, the potential of each In atom in 
the structure is obtained as: 

ji ^ — fi 
Win (As -u In (InAs) H — f; In (InP) (A4) 
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